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We consider an ultracold rotating Bose gas in a harmonic trap close to the critical angular ve- 
locity so that the system can be considered to be confined to the lowest Landau level. With this 
assumption we prove that the Gross-Pitaevskii energy functional accurately describes the ground 
state energy of the corresponding A''-body Hamiltonian with contact interaction provided the total 
angular momentum L is much less than A^^. While the Gross-Pitaevskii energy is always an obvious 
variational upper bound to the ground state energy, a more refined analysis is needed to establish 
it as an exact lower bound. We also discuss the question of Bose-Einstein condensation in the pa- 
rameter range considered. Coherent states together with inequalities in spaces of analytic functions 
are the main technical tools. 



I. INTRODUCTION 

A Bose gas rotating in a harmonic trap has a criti- 
cal rotation speed above which the trap cannot confine 
it against centrifugal forces. If the trapping potential 
equals ^{oj^r"^ -\- lo^x\), with m the particle mass and 
r = ^Jx\+ x\ the distance from the axis of rotation, 
then the critical angular velocity is cjj^. In a reference 
frame rotating with angular velocity f2 the Hamiltonian 
for one particle is 



(1) 



with L the component of the angular momentum along 
the rotation axis. It is convenient and instructive to com- 
plete the square and write ^ as 

±.ip-Af + '^xl + {^^-n)L (2) 

where A = mLUj_{x2, —xi, 0). In the rapidly rotating case, 
where < U!± — Q <C minjtjj^, lu^}, it is natural to restrict 
the allowed wave functions to the ground state space of 
the first two terms in ([2]), which we denote by H, and 
this restriction will be made in this paper. The space Ti. 
consists of functions in the lowest Landau level (LLL) for 
motion in the plane perpendicular to the axis of rotation, 
multiplied by a fixed Gaussian in the xs-direction. Apart 
from the irrelevant additive constants 2huj± and huj^, the 
kinetic energy in Ti is simply 



LUL 



(3) 
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with Lu = uj± — n > 0. Note that L is non-negative for 
fimctions in Ti.. Its eigenvalues are 0, h, 2h, .... 

To characterize the functions in the space Ti., it is nat- 
ural to introduce complex notation, z ^ xi + 1x2- Func- 
tions in TL are of the form 



/(^)exp 



2h 



muj3 2 



(4) 



with / an analytic function. All the freedom is in / since 
the Gaussian is fixed. If the trapping potential in the X3- 
direction is not quadratic, the Gaussian in the x^-variahle 
has to be replaced by the appropriate ground state wave 
function. 

A fancy way of saying this is that our Hilbert space 
Ti consists of analytic functions on the complex plane C 
with inner product given by 



0(z)Xz)e-l^l' dz, 



(5) 



where dz is short for dxi dx2- For simplicity we choose 
units such that m — h — = 1. The eigenfunction 
of the angular momentum L corresponding to the eigen- 
value n is simply z". In other words, L = zdz- We 



remark that the expectation value of 



P in this 



state is n + 1. 

For a system of N bosons, the appropriate wave func- 
tions are analytic and symmetric functions of the bosons 
coordinates zi, . . . , zn. The Hilbert space is thus = 
®symm^- The kinetic energy is simply uj times the total 
angular momentum. 

In addition to the kinetic energy, there is pairwise in- 
teraction among the bosons. It is assumed to be short 
range compared to any other characteristic length in the 
system, and can be modeled by gJ2i<i<j<N ^i-^i ~ ^3) 
for some coupling constant g > 0. Physically, this g is 

1 /2 

proportional to awg where a is the scattering length of 
the three-dimensional interaction potential. On the full. 



2 



original, Hilbert space (^^^^L"^ (M.'^) , a 5-function as a 
repulsive interaction potential is meaningless. On the 
subspace Tijv the matrix elements of S{zi — Zj) make per- 
fect sense, however, and define a bounded operator Stj. 
Using the analyticity of the wave functions this operator 
is easily shown to act as 

{Sl2^) {Zi, Z2) - 2^V'((^1 + 22)/2, (Zl + Z2)/2) , (6) 

which takes analytic functions into analytic functions. Its 
matrix elements in a two-particle function 'ip{zi, Z2) are 



(7) 



The dimensional reduction from three to two dimen- 
sions for the A^-body problem and the restriction to the 
LLL is, of course, only reasonable if the interaction en- 
ergy per particle is much less than the energy gap 2hujj_ 
between Landau levels and the gap huj3 for the motion 
in the cca-direction. For a dilute gas the interaction en- 
ergy per particle is of the order ap where p is the av- 
erage three-dimensional density Provided Ng/uj is 
not small we can estimate p by noting that the effec- 
tive radius R of the system can be obtained by equating 
ap ~ Na/ {R^oj^^^'^) ~ Ng/R^ with the kinetic energy 
LoL ~ LoR^. This gives R ~ {Ng/ajY^^ and the condition 
for the restriction to the LLL becomes 



mm{uj±,uJ3}. 



(8) 



The physics of rapidly rotating ultracold Bose gases 
close to the LLL regime has been the subject of many the- 
oretical and experimental investigations in recent years, 
starting with the papers 0, 0, U H, @] ■ The recent re- 
views 0, H, [13] contain extensive lists of references on 
this subject. On the experimental side we mention in 
particular the papers [llf, [l^, [l3j] that report on exper- 
iments with rotational frequencies exceeding 99% of the 
trap frequency. 



II. MODEL AND MAIN RESULT 

The discussion in the Introduction leads us to the fol- 
lowing well-known model (see, e.g., 0, 0, for ^ 
bosons with repulsive short-range pairwise interactions: 



N 



9 

l<i<j<N 



(9) 



It acts on analytic and symmetric functions of N vari- 
ables S C. The angular momentum operators are 
Li ~ Zidzi, and Sij acts as in (jS]). The parameters uj 
and g are assumed to be positive. The rigorous deriva- 
tion of this model from the 3D Schrodinger equation for 
particles interacting with short, but finite range poten- 
tials will be presented elsewhere (141] . 



li<j 

course, 
five. 



Note that the two terms in H commute with each 
other, and hence can be diagonalized simultaneously. 
An exact eigenfunction with total angular momentum 
Ltot = J2f=i = N{N — 1) is the Laughhn state 
J|j^^-(zi — ZjY for which the interaction vanishes, of 
For smaller L the interaction is strictly posi- 
It is expected that the lowest eigenvalue of the 
interaction is of the order N when Ltot is of the order 
N"^ . The lower boundary of the joint spectrum of inter- 
action and angular momentum has been called the yrast 
curve [15], a term originating in nuclear physics. Upper 
bounds on this curve can be found by variational calcula- 
tions, and a host of trial functions with many interesting 
properties related to the Fractional Quantum Hall Ef- 
fect (FQHE) have been employed for this purpose, see, 
e.g., B, M [il [11 S HH. A challenging problem is to 
establish the missing reliable lower bounds. 

In this paper we consider the case when Ltot is much 
smaller than iV^, i.e, we investigate one corner of the 
asymptotics of the yrast region. In this case, we shall 
show that the ground state energy of H is well approxi- 
mated (rigorously so in the N ^ 00 limit) by the Gross- 
Pitaevskii energy. The latter energy is the minimum of 
the Gross-Pitaevskii functional 



£GP(0) = c.(0|L|0) + ^ / |0(z)|4 e-2|-l' dz (10) 



9 
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over all analytic functions with (0|0) — N. This can 
also be viewed as a Hartree approximation to the many- 
body system, where one takes the expectation value of H 
with a simple product function 11^=1 4'{zi) (and ignores 
a factor {N — 1)/N). The minimization problem for (fTU]) 
and the properties of the minimizers, in particular their 
vortex structure, have been studied in some detail by 
many authors including [1, [H [H HI HI ■ 

To state our result precisely, let EQ{N,uj,g) be the 
ground state energy of H, and let E'^^{N,uj,g) be the 
Gross-Pitaevskii energy. We will find a positive, finite 
constant C such that 



E'^''iN,u;,g)>EoiN,u;,g) 
> E'~^^{N,Lj,g) 



1-cm 

\NujJ 



1/10' 



(11) 



for all TV, uj and g such that gN/ui is bounded below 
by some (arbitrary) fixed constant. More precisely, for 
any c > there exists a C < 00 such that (|lip holds if 
gN /uj > c. In particular, this implies that the ratio of 
Ea and L'^^ is close to one if g <C Nuj and g > N^^uj. 

Note that, by simple scaling, E'^^ {N,uj, g) — 
NujE'-^^ {1,1, gN/uj). Hence the contrary case of small 
gN/uj is not particularly interesting; as gN/uj — s- one 
obtains a non-interacting gas. On the other hand, if 
gN /uj is not small the kinetic and interaction energy 
are of the same order of magnitude, which, according 
to the previous back-of-the-envelope analysis, implies 
E'^^ {N,uj, g) ^ Ny/Ngui. Hence the total angular mo- 
mentum, which is obtained by taking the derivative of 



3 



the energy with respect to w, is of the order N Ng/u, 
which is much less than if and only if g ^ Nui. 
Hence our bound (jlip covers the whole parameter regime 
itot < N\ 

We note also that in terms of the filling factor v ~ 
N^/{2Ltot) 6] the parameter regime g/{Nuj) <C 1 corre- 
sponds to 3> 1. In contrast, the Laughlin wave function 
has filling factor v = 1/2 for TV large. Between these two 
extremes rich physics related to the FQHE is expected 
[lol | , but this regime is apparently still out of experimen- 
tal reach. 

Before giving the proof of pT|) . we shall discuss some 
of its implications for the yrast line. Recall that the 
yrast energy Io{L) is defined as the ground state energy 
of X)i<j ^ij sector of total angular momentum L. 

What pT|) says is that 



LuL + gIo{L)>E''''{N,iu,g) 



\NujJ 



1/10' 



(12) 



for any i > 0. As Ng/uj gets large, it is expected that 



E^^{N,g,uj) « bN^Nguj 



(13) 



to leading order, with h ~ 0.57. An upper bound with this 
value of b was actually derived in [26i| , but a lower bound 
with the same b is still open. We remark that it is easy 
to derive a lower bound with b = -^/S/Ott w 0.53. This 
follows using = / |0(z)|2(|z|2 - l)e-l^l'dz and 

minimizing /(u;|z|2|</,(z)|2e-l^l' +(.g/2)|0(z)|4 e'^l^l' )dz 
over all dropping the LLL condition. Thus, 0.53 < 
b < 0.57. Using our bound we deduce from (fT^ that 



loiL) 



> 



b^N^ 
4 L 



for TV < L < iV^ 



For 2 < L < iV it is well known that Io{L) = 
N{At:)-^{N - 1 - L/2) [3, i, 113, [H, m. Our bound 
(fT5)) reproduces this result for large N since, as we shall 
prove below, 



Ifg 
in 



if 



< 



(14) 



which is just the GP energy of the constant function. In 
particular, the zero angular momentum state is a mini- 
mizer of the GP functional for gN < Sttuj. Eqs. (fT^ and 
(fTi)) imply that 



loiL) > 



N 
An 



N 



1-C 



8n 



1/10' 



(15) 



for all L > 0. In order to prove (|T4l) . we note that the GP 
energy (jlOp equals a certain two-particle energy, namely 

^°^('/') = \{<t>®(t>\^{Li + L2)+g5i2\(l)®(l)) ■ 



The operator 5i2 commutes with Li + L2, and satisfies 
6I2 = (27r)~^5i2. That is, 2n6i2 is a projection operator. 
We claim that 



> 4. 



1 



2nSi2\ (j) (/)) 
(16) 

for any (j> £ H. This clearly implies (fT4|l . To see (fT6|) . we 
can decompose the two-particle function 0(zi)0(z2) into 
a sum X]n>o ("^1 ' •^2) of functions of given total angular 
momentum n. For the terms with n > 4 the bound (jl6p 
certainly holds. Also for n = and n = 1 it holds, since 
1 — 27r5i2 = in this case. We thus have to consider 
only the cases n — 2 and n = 3. On the subspace with 
n = 2, (1 — 2n5i2)4i® (p is proportional to the normalized 
function 772(2:1, Z2) = (zi — Z2)^/\/8- It is straightforward 
to check that \{4> ® 4>\ri2)\ < ||''/'2||/'\/2 which implies the 
desired result. Similarly, on the subspace with n = 3, (1— 
27r(5i2)0(8)</' is proportional to 773(21, Z2) = (zi — 22)^(^1 -t- 
Z2)/VlQ- One checks that \{cj) <^ < IIV'sllVsTi, 

which is what is needed. 

Everything now depends on the proof of (fTTj) and the 
most of the rest of this paper is devoted to this task. 
The first inequality in (fTTj) follows easily by taking prod- 
uct wave functions as trial wave functions, as mentioned 
earlier. Hence it remains to show the second inequality 
in (fTTj) . which is the lower bound on i?o. 

For the proof we shall employ the technique of coher- 
ent states and c-number substitutions 27, 28] that was 
used in [29| to solve a related problem, namely to derive 
the Gross-Pitaevskii equation from the three-dimensional 
many-body Hamiltonian of a rotating Bose gas away from 
the critical rotational frequency. This was done in the 
Gross Pitaevskii limit where both w and Ng are fixed 
and order unity. The situation discussed in the present 
paper is partly simpler than that in [2£|] , where the inter- 
action was described by an arbitrary repulsive potential 
of short range instead of a contact interaction. On the 
other hand we now have to face new problems. In the 
present setting the coupling parameter gN/co can vary 
with iV (as long as it is <C N'^) and this fact requires 
considerably more delicate estimates on matrix elements 
of the interaction between two-particle states than in [2^ . 
In particular, our bound is asymptotically exact in the 
whole "Thomas-Fermi" regime. In order to make the 
proof more transparent it will be divided into seven steps. 



III. DERIVATION OF THE MAIN INEQUALITY 

m 

A. Step 1 

As a first step towards deriving a lower bound on the 
ground state energy of H, we consider a slightly big- 
ger Hamiltonian H' which is constructed in the following 
way. The reason for introducing H' is that the addi- 
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tional small positive term is needed to compensate neg- 
ative terms that will occur at a later stage of the proof. 
For integers J, let 



xAL) 




L-y/jTT) 



for L<J 
for L > J + 1. 

(17) 

Pick an integer Jq > 1. We first show that there exists 
a J with Jo < J < 2 Jo such that 



N 



(18) 



where ( • ) denote the expectation value in the ground 
state of H. To see this, note that 



J2xjiL) < 1 + 2 / exp 



J>0 



s ds 



< 1 + 4V7ri . 



In particular, since xj(0) — 0; we have J2.jXj{L) < 
(1 + 4y^)-\/L for all integers L > 0. This implies that 



2Jo-l 



^ E XjiL) < (1 + 4V^)^ < 2^/2(1 + 4V^) 



Jo 



J3/2 ' 



where we have used that L > Jq in order that the left 
side be non-zero, and Jq > J/2. The expectation value 
of is bounded by the total energy divided by uj. 

Hence the bound (fTS]) holds on average for Jo < J < 2 Jq 
and must thus hold for at least one such J. 

We can now pick an 77 > and consider the modified 
Hamiltonian 



N 



i=l 



By (jlSp its ground state energy is bounded from above 
by £^0(1 + br]), with b = 2^/2(1 + 4^). In other words, 
Eq is bounded from below by the ground state energy of 
H' divided by (1 + brj). We choose 77 = {g/Nujy^^^. 



B. Step 2 

In order to derive a lower bound on H' , we shall first 
extend it to Fock space in the usual way. We do this in 
order to be able to utilize coherent states and the lower 
and upper symbols of the extended operator. 

Let fj{z) = z-' (ttj!)"^/^ denote the normalized eigen- 
functions of the angular momentum operator L with 
eigenvalue j G N. Let al and denote the correspond- 
ing creation and annihilation operators. On Fock space. 



consider the operator 



i>0 



ijk 



(19) 



for some /i > 0. For simplicity we denote 612 simply by S. 
Note that because of conservation of angular momentum 
all other matrix elements of S vanish. The Hamiltonian 
H' can be viewed as the restriction of H to the sub- 
space containing exactly N particles. The value of /i is 
irrelevant, since the term multiplying fi vanishes in this 
subspace. In particular, the ground state energy of H' 
is bounded from below by the ground state energy of H, 
for any value of fi. 

We introduce coherent states for all angular momen- 
tum states up to J. That is, in Eq. (|19[) we normal-order 
all the a| and in the usual way (which is only rele- 
vant for the term multiplying /i since all other terms are 
already normal ordered) and then replace all by com- 
plex numbers Q and all a^- by the conjugate numbers Q 
for < i < J. The resulting operator on the Fock space 
generated by the modes > J is called the lower sym- 
bol of H and will be denoted by hi{(), where ( stands 

for (Co, ■ • • ,0)- Note that hg{() does not conserve the 
particle number, which explains why it was necessary to 
embed our iV-particle Hilbert space into Fock space. 

The ground state energy of H is bounded from below by 
the ground state energy of the upper symbol hu{C) when 
minimized over all parameters C- The upper symbol is 
given in terms of the lower symbol as [l^l 

(20) 

where Ui{C) - Z^i^o ^'^C* ^'^(C) equals 
-oj'^i-2g ^ {(l>c <E) (pi\S\(l)i; (g) (fi^) 

i<J 0<i<.J 



0<i<,J k>J 

fi [(2||0<;||2 -2N + 2N> + 1) (J + 1) + 2\\(l)^f] 



and 



hdC) 



0<i,j<J 

Here, 0^ = E/=oO'/'.. UcW = (^cl-^c) = Eto iGf 
and = Efc>j'^]'^j- Note that in a Taylor expan- 
sion of the exponential in ((20|) only the first three terms 
contribute since /if (C) is a polynomial of degree four. 
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For a lower bound, we can use the fact that U2{C) > 0. 
Moreover, to bound the various terms in Ui{C) we can 
use the fact that X]i<j ^ Tr^^e'^l". It is then 

easy to see that Ui {() is bounded from below as 



U! 



u,{C)>~-JiJ- 

-2N' 



1)-('^ + 2m(J 



2) m\ 



1) + ^ 

IT 



(21) 



In order for the first term to be much smaller than the 
GP energy E'^^{N,g,uj) ~ N^/Ngcu it is clear that we 
can not replace all modes by a c-number but must rather 
require that J <^ N{g/NLjy/^. 



C. Step 3 

As we have explained in the previous step, H' is 
bounded from below by the minimum over C of the 
ground state energy of huiC)- Suppose the minimum is 
attained at some Co: and let ( • )o denote the correspond- 
ing ground state expectation value in the Fock space of 
the modes > J. The next step is to use a simple lower 
bound on /lu(Co) to bound the 77-dependent term in H 
involving the modes > J. 

Let us choose the parameters fi and J in following 
way. Recall from step 1 that Jq < J < 2Jo, with Jo 
an arbitrary integer that we are free to choose. We take 
Jo = [N{g/NLj)^^^^\, where [t\ denotes the integer part 
of a number t > 0, and /i — loJ/{AN). The propor- 
tionality constants here are chosen more for convenience 
than for optimality. Note that Jq is a big number if 
g/{Nu}) < 1 and gN/uj > 1. 

Recall that rj = {g/Nivy/^^. We claim that, for small 
g/Nu;, 



In particular, using (|2T|) . 



i>J 



25 



(J + 1){lu - 2//) - — - 2/iiV 

TT 



-N['^ + 2^(J + 2)) - ^J{J + 1) 



^ + 2^,iJ + 2)]iU^r-N). 

TT 



With the aid of Schwarz's inequality, we obtain 



i>J 



N 



25 



{J + l){uj - 2ii) - — - 2nN 

TT 



2g 



2^(J + 2) 



uj , , (^-|-2^(J + 2))' , , 
--J(J+1)-1^ ^. (24) 

For our choice of /i and J, /i <C 1^ and g ojJ for 
small g/Ncu. Moreover, iiN — uoJ/A. Hence the second 
term on the right side of ([M)) is positive for small g/Noj. 
Moreover, since E*^^ {N,uj, g) ^ Ny/Nujg for Ng > w, 
as remarked earlier, all the terms in the last line of ([^^ 
are much smaller than E'~^^ {N , lu , g) if 5 ^ Nu! and, in 
particular, are bounded from below by —E'^^/2. Since 
the ground state energy of hu{Co) is bounded above by 
Eo{N,uj,g)il + brj) < 1 .5 E'=^^ {N , lj , g) , Inequality ^ 
follows. 



\i>J I 



E'^^{N,uj,g) 



(22) 



This bound is similar to (|18p except for the additional 
factor 77 on the left side and E'^^ replacing Eq on the right 
side. In order to prove (P^ . we can use the positivity of 
the interaction to obtain the following lower bound on 
/i,j(C). The lower symbol h^^Q) is bounded from below as 

hfXO>vY.^ji^)a\a,+Lo{J+l)N> 



i>.j 



{\\^^\\^-Nf + 2N>{U^r-N) 



+ {N>r+\m 



(23) 



D. step 4 



We now give a more refined lower bound on he((^). In- 
stead of dropping all the interaction terms, as we did in 
the previous step in order to obtain the bound (j23p . we 
use the fact that 



S>Pj(E) PjSPj (E)Pj + Pj (E) Pj5{Pj ® Qj + Qj (E) Pj) 
+ (Pj (S>Qj + Qj(S) Pj)5Pj (g) Pj , (25) 

where Q,j = 1—Pj and Pj denotes the projection onto the 
J -|- 1 dimensional subspace of 7i spanned by eigenfunc- 
tions of the angular momentum with L < J. Inequality 
(|25p is a simple Schwarz inequality and follows from posi- 
tivity of S and the fact that Pj®Pj5Qj®Q,j = 0. Using 
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(|25)) . we see that /i^(C) is bounded from below by 

j>J 

{Ucr-Nf + 2N>{Ucr-N) 
+ 3 ^((/>C ® 4'c\S\4'C ® Vk)ak + h.c. . 

k>J 

For the remaining terms in the upper symbol /iu(C) 
proceed as in the previous step. We retain — 
NY /2 for later use, however, and arrive at 

k>J 

-n(^ + J + 2) j - I J( J + 1) 



(f + 2MJ + 2))^ 
2fi 



(26) 



The right side of (pS)) contains the desired quantity 
f ^^(^^). The second term guarantees that ||'/>ij |P is close 
to N. All the terms in the last two lines are small com- 
pared to E'^^{N,uj,g) for our choice of /x and J. Hence 
we are left with giving a bound on the third term in ([26)) . 
which is linear in creation and annihilation operators. It 
is this bound that requires most effort; it will be com- 
pleted in the remaining three steps. 



E. Step 5 

A simple Schwarz inequality shows that, for any se- 
quence of complex numbers Ck and positive numbers e^, 



E {^^-k + 44) < E ^ + E ^^4«. ■ (27) 

We apply this to our case, where Ck = — (0c ® '/'cl'^l'^C ® 
(fik) for k > J, and = otherwise. We pick some k > 
and choose ej, = K^^xj(fc) for k > J, with Xj(^) given 
in ini). Then 



E^ = «('/'ckx7'(^KI'/'c>, (28) 
k ^'^ 

where denotes the multiplication operator P(^{z) = 
|(/)c(2)pe-l^l' and 

OO ^ 

x:j^{l)= E — 7^Tl'^^)(^fcl• 



(For simplicity, we abuse the notation slightly, since 
does not leave 7i invariant; {'^k\pc,\4'c) makes perfectly 
good sense, however.) 

In Step 7 below we shall show that 



3/2 



(29) 

This inequality may seem surprising at first sight. The 
function P(i\4>c) contains angular momenta up to L < 2 J, 
and xJ^(2J) grows exponentially with J. As will be- 
come clear below, however, |(</'fe|Pcl'?^c)P exponentially 
small unless k is smaller than roug hly J ^ y/j- bmce 
X]^{J + ^fj) is bounded independently of J, this makes 
plausible. 
Ahogether, dUD-dlHl) imply that 



3/2 



E ^'=4 



E (^'^"fc + 44) 

k 

<6k / |0c(z)|4e-2|^l'dz 
We have shown in Step 3 above that 

\E^^°^"7/' ,^,j3/2 

in the ground state of /im(Co)- After optimizing over k 
we get the following lower bound on — X]/c('^'='^fe + 4'^!) 
(valid as an expectation value in the ground state): 



-4V3 



|0c(^)|''e-2|^l'dz 



Finally, using Holder's inequality this expression is 
bounded from below by 



for arbitrary [3 > Q. 



F. Step 6 



The analysis in the preceding steps has shown that 

K{Qo) > (1 - m^'^'iK) + f {WKf - Nf 

_ iV Ji + 2p{J + 2) - ^ 

V TT / 2p 



65 /^GP(7V^y 



(30) 
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As mentioned above, we choose /i = LoJj (4iV) and 
J - N{g/Nujf/^^. Moreover, we make the choice 
rj = (3 — {g/NtoY^^^. Then all the error terms are small 
compared to E'^^ {N , uj , g) , namely they are of the order 
(5/iVcc;)i/i0i;GP. 

Since the ground state energy of /i(Co) is bounded 
above by E'^^{N,u;,g){l + hr^) < l.b E^^ {N,uj, g) and 
all the error terms in ([50]) are bounded by E'^^ /2 for 
small g/{Nuj), we see that 



\<t>c\\ 



N < 



2E^^{N,Lu,g) 



1/2 



\NloJ 



1/10 



for C = Co a-iid hence 



1/10 



\NljJ 



for some constant C > 0. This yields the lower bound in 



G. Step 7 

The only thing left to do is to prove Inequality ([29l 
For convenience, we shall introduce the notation 



|0(z)|Pe-(f/2)l^l'dz 



i/p 



for 1 < p < oo. We can bound 

{^c\pcXj\l)pM < Uc\\1\\VPCXj\L)^\\ . 

The latter norm is an operator norm on L^{C, e~'^' dz) 
(defined as the maximum expectation value), and equals 

the operator norm of Xj^ {L)PC\J Xj^iL) on TC. 

Next, we derive a pointwise bound on P({z) = 

note the integral kernel of the projection onto the sub- 
space of H with angular momentum L < J. With the 
aid of Holder's inequality. 



The LP norms with p ^ 2 are not given by simple expres- 
sions, however. In the following, we will show that they 
can be bounded above by a quantity independent of J. 
More precisely, for any 1 < p < oo, we shall show that 



\\PAz,-)\\p< 



i/p 



TT sin(7r/p) 



independently of J. In order to see this, note that, for 
fixed z and \z'\, the function 9 i-^ Pj{z,\z'\e^^ ) is ob- 
tained from Poo (2, \ z'\ e*^ ) by restricting the Fourier com- 
ponents to < j < J. This restriction is a bounded 
operation (uniformly in J) on LP{T) for 1 < p < 00 j30l . 
Ch. II, Sec. 1], and hence \\Pj{z, ■ )\\p is bounded by a 
constant times ||Poo(2, Ollp- fact, an upper bound 
on the optimal constant in this inequality is given by the 
norm of the Riesz projection (viewed as an operator from 
LP{T) to itself), which is known to equal 1/ sin(7r/p) [3]| . 
Hence we have the bound 



\\PA^, ■)\\p 



< 



sin(7r/p) 
i2n/py/P 



Poo (2 J ' 

\zf/2 



TT sm(n/p) 



(31) 



The last equality follows from the fact that Poo(z, z') 



Let Cp denote the constant in square brackets in (j3ip 
taken to the power 2p/(4 — 2p). We have thus shown 
that, for 1 < p < 4/3 and a{p) = (4 - 3p)/(4 - 2p) 



(^)<Cp||0clll(E,=ol^^(^)l'""'"' 



i(p) 



A possible choice is p = 6/5, in which case Cp = 
(10/3)V2(5/37r)i/4 ~ 4.02 and a{p) = 1/4. Hence, for 
any function \ip) = J2k>j dkl'fik), we have that 



\/x7WpcVx7(I) 



i^f ^'lYl ^ <4.O2||0cllI E |4pexp[i(x/I-x/jTT) 



k>J+l 



Pj(z,z')4>dz')e-\''\'' dz' 



zll^e-l^l^ 



1/4 



e-l^l dz. 
(32) 



< ||0cll4||P7(^, OIU/S. 

Another application of Holder's inequality yields 



Pj{z, 



1 4/3 



< \\Pj{-^ 



j(4-3p)/(4-2p)| 



Pj{z, 



|p/(4-2p) 



for any 1 < p < 4/3. 

The LF' norm oi Pj{z, • ) is given by 



Here we have used the fact that our upper bound on P(^{z) 
is radial, i.e., depends only on |z|, and hence there are 
no off-diagonal terms on the right side. An application 
of Jensen's inequality for the concave function t t-^ t^/^ 
implies that 



\pjiz, 



< 



1/4 
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From Stirling's formula [S^l it follows that j\ > 
{j/ey^yi + j for any j > 0. Using this bound it is easy 
to see that 



|^,(z)|2e-l^l^ < 



1 



Similarly, 

7r\/8e \z\ 

for fc > 1. Hence 



< 



33/2 



1 



TrV2e VI + j Jo 



dx 



< 



33/2 



1 



/47re a/I + j 



To obtain the last inequality, we simply extended the 
integral to the whole real axis. The sum over j can be 
bounded from above by the integral 



-(V7-Vfe)'/2 



3=0 

Jo Vs 



< 2 

< 2 



-00 
00 



-{t+^/J+l-Vkf/2 ^-~t{Vk-^/J+l) 



Stt e- 



-(Vfc-Vj+T)V2 



The one-fourth root of this expression cancels exactly the 
exponential factor in p2p . and we arrive at 



the GP functional pO|) due to breaking of rotational sym- 
metry makes the precise statement a little complicated, 
but in brief the result is as follows: If ^'jv is a sequence of 
ground states of ([5]) for iV = 1, 2, . . . and are the cor- 
responding normalized 1-particle density matrices, then 
any limit 7 of the sequence as — > 00 with the 
coupling constant gN /lij fixed is a convex combination of 
projectors onto normalized minimizers of the GP func- 
tional (fTO| for the given coupling constant. A precise 
formulation is given in Theorem 2 in (29j . A key step in 
the proof is to extend the bounds on the ground state 
energy of 7? to perturbed operators H^^'> = H + J^iLi Si 
where 5 is a bounded operator on the 1-particle space Ti. 
and Si the corresponding operator on the i-th factor in 



^syinm^- Thcsc bounds lead to the inequalities 



< TrS'7 < max(0|S'|0) 



(33) 



for any limit 7 of the sequence of density matrices jn, 
where the max resp. min is taken over all normalized min- 
imizers 4> of (|10p with coupling parameter gN/u. With 
the aid of some arguments from convex analysis one can 
then conclude as in 29] that 7 has a representation in 
terms of projectors onto normalized minimizers (p 

of the GP functional. The appropriate extension of this 
result to the Thomas-Fermi limit is still an open problem. 



V. CONCLUSION 

We analyzed the ground state energy of a rapidly rotat- 
ing Bose gas in a harmonic trap, under the usual assump- 
tion that the particles are restricted to the lowest Landau 
level, that the two-body interaction is a (5-function, and 
that the number of particles, TV, is very large. For the 
limiting situation we prove, rigorously, that the Gross- 
Pitaevskii energy, appropriately modified to the lowest 
Landau level, is exact provided g ^ Nui. Here, g is the 
strength of the (5-function interaction, and uj is the dif- 
ference between the trapping and rotation frequencies. 



IV. BOSE-EINSTEIN CONDENSATION 

The technique employed for the energy bounds can be 
used to show that the model ^ exhibits Bose-Einstein 
condensation in the ground state if — > cxd with gN/uj 
fixed, in complete analogy to the proof of the correspond- 
ing result in ^2^] . The nonuniqueness of the minimizers of 
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